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Abstract. Lya hlobs (LABs) are spatially extended Lya nebulae seen at high redshift. The origin of Lya emission in the LABs 
is still unclear and under dehate. To study their heating mechanism(s), we present Australia Telescope Compact Array (ATCA) 
ohservations of the 20 cm radio emission and Herschel PACS and SPIRE measurements of the far-infrared (PIR) emission 
towards the four LABs in the protocluster J2143-4423 at z=2.38. Among the four LABs, B6 and B7 are detected in the radio 
with fluxes of 67±17 /rjy and 77±16 /rJy, respectively, and B5 is marginally detected at 3 cr (51±16 yuJy). For all detected 
sources, their radio positions are consistent with the central positions of the LABs. B6 and B7 are obviously also detected in 
the FIR. By fitting the data with different templates, we obtained redshifts of 2.20)jj35 for B6 and 2.20)|5'3 q for B7 which are 
consistent with the redshift of the Lya emission within uncertainties, indicating that both FIR sources are likely associated with 
the LABs. The associated FIR emission in B6 and B7 and high star formation rates strongly favor star formation in galaxies 
as an important powering source for the Lya emission in both LABs. However, the other two, B1 and B5, are predominantly 
driven by the active galactic nuclei or other sources of energy still to be specified, but not mainly by star formation. In general, 
the LABs are powered by quite diverse sources of energy. 
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1. Introduction 

High-redshift star-forming galaxies are becoming an impor¬ 
tant probe of galaxy formation, reionization and cosmology 
(Robertson et al. 2010; Shapley 2011). A popular method for 
finding high redshift star forming galaxies is to target their of¬ 
ten bright Lya emission (Partridge & Peebles 1967). This emis¬ 
sion can be easily detected in narrow-band imaging surveys, 
and can be further confirmed by spectroscopic observations 
(Hu et al. 1998; Ouchi et al. 2008; Yamada et al. 2012a,b). In 


* Herschel (Pilbratt et al. 2010) is an ESA space observatory with 
science instruments provided by European-led Principal Investigator 
consortia and with important participation from NASA. 

** email: yiping.ao@nao.ac.jp 


addition to discovering numerous Lya emitters (LAEs), a par¬ 
ticular class of objects, also known as ”Lya blobs” (LABs), has 
been most commonly found in the dense environment of star¬ 
forming galaxies at high redshift, and they are very extended 
(30 to 200 kpc) and Lya-luminous (10"^^ to 10“*"^ erg s“^) (see, 
e.g., Francis et al. 1996; Steidel et al. 2000; Palunas et al. 2004; 
Matsuda et al. 2004, 2009, 2011; Dey et al. 2005; Saito et al. 
2006; Yang et al. 2009, 2010; Erb et al. 2011; Prescott et al. 
2012a, 2013; Bridge et al. 2013). In contrast to the large Lya 
nebulae surrounding some high-redshift radio galaxies (e.g., 
Reuland et al. 2003; Venemans et al. 2007), these objects do not 
always have obvious sources for energy responsible for their 
strong emission. 
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While the LABs’ preferential location in overdense en¬ 
vironments indicates an association with massive galaxy for¬ 
mation, the origin of Lya emission in the LABs is still un¬ 
clear and under debate (Faucher-Giguere et al. 2010; Cen & 
Zheng 2013; Yajima et al. 2013). Proposed sources have gen¬ 
erally fallen into two categories: cooling radiation from cold 
streams of gas accreting onto galaxies (e.g., Haiman et al. 
2000; Dijkstra & Loeb 2009; Goerdt et al. 2010) and photoion¬ 
ization/recombination from starbursts or active galactic nuclei 
(AGNs) (e.g., Taniguchi & Shioya 2000; Furlanetto et al. 2005; 
Mori & Umemura 2006; Zheng et al. 2011). Supporting evi¬ 
dence for the cooling flow scenario comes from those LABs 
lacking any visible power source (e.g., Nilsson et al. 2006; 
Smith & Jarvis 2007). Ionizing photons from young stars in 
star-forming galaxies and/or AGNs can ionize neutral hydrogen 
atoms and the subsequent recombination gives off Lycr emis¬ 
sion. The resonant scattering of Lya photons in the circum- 
galactic medium makes the emission extended (Geach et al. 
2005, 2009; Colbert et al. 2006, 2011; Beelen et al. 2008; Webb 
et al. 2009; Zheng et al. 2011; Cen & Zheng 2013; Overzier et 
al. 2013). 

Except for cooling flows and photoionization from star¬ 
forming galaxies and/or AGNs, other possible mechanisms, 
such as galactic super-winds and obscured AGNs, are also pro¬ 
posed to explain the nature of the LABs (e.g., Ohyama et al. 
2003; Wilman et al. 2005; Colbert et al. 2006; Matsuda et al. 
2007). All these sources of energy may be activated in an en¬ 
vironment where violent interactions are frequent between gas 
rich galaxies as expected in over-dense regions at high redshift 
(Matsuda et al. 2009, 2011; Prescott et al. 2012b; Kubo et al. 
2013). 

The 110 Mpc filament with 37 LAEs related to the pro¬ 
tocluster J2143-4423 at z=2.38 (Erancis et al. 1996, 2004; 
Palunas et al. 2004) is one of the largest known structures 
at high redshift, and this field also includes four large ex¬ 
tended LABs with extensions of ~ 50 kpc and above, named 
Bl, B5, B6 and B7. In this paper, we present our deep ra¬ 
dio observations and Herschel released far-infrared (EIR) data 
in J2143-4423 to study the powering source of these LABs. 
Throughout this paper, we use a A cosmology with Ho = 
67.3 km s ' Mpc ', Ga = 0.685 and = 0.315 (Planck 
Collaboration XVI 2013), and 1" corresponds to 8.37 kpc at 
z=2.38. 

2. Observations 

2.1. ATCA observations 

We observed J2143-4423 with the Australia Telescope 
Compact Array (ATCA|^ in its extended configuration 6A. 
During the observations from 2009 June 14 to 17, only five 
out of six antennas were available. The observations were 
performed at a central frequency of 1.75 GHz. We used the 
Compact Array Broadband Backend (Wilson et al. 2011) in a 

' The Australia Telescope Compact Array is part of the Australia 
Telescope, which is funded by the Commonwealth of Australia for 
operation as a National Facility managed by CSIRO. 


wide-hand mode, with a total bandwidth of 2 GHz and a chan¬ 
nel width of 1 MHz. The nearby source PKS 2134-470 served 
as a gain calibrator. Absolute fluxes were calibrated with the 
ATCA standard PKS 1934-638. The total observing time was 
about 70 hours. 

The data were reduced with the MIRIAD software pack¬ 
age. Although the observations were carried out with a to¬ 
tal bandwidth of 2 GHz, the effective bandwidth was about 
489 MHz with a central frequency of 1.51 GHz. We carefully 
flagged the channels affected by radio frequency interference 
(REI) by checking the visibility data sorted by time, channels 
and baselines. The image was deconvolved with MIRIAD task 
MECLEAN, and task SELECAL was used to reduce the noise 
from strong radio continuum sources. We first created cleaned 
images in a normal procedure and made model images for 
the strong sources. The models were used as inputs for task 
SELECAL to perform self-calibration of visibility data. We ran 
this cycle for three times, and then obtained the model images 
to create the visibility data with self-calibration, which were 
used to make the hnal images. The noise of the images after ap¬ 
plying self-calibration was about one order of magnitude lower 
than that without self-calibration. The field of view was about 
31 arcmins and the synthesized beam size was 7.8"x4.8". The 
noise was about 15 juJy/beam before applying primary beam 
correction. 

2.2. Archival Herschel observations 

Herschel observations towards J2143-4423 were carried out 
with PACS (Poglitsch et al. 2010) at 100 and 160 jum and 
SPIRE (Griffin et al. 2010) at 250, 350 and 500 yum in 2010 
to 2011. J2143-4423 was imaged in a held size of 15'xl5' for 
each band, and the observing time was ~2.9 hours for PACS 
{Herschel OD; 686) and ~0.6 hours for SPIRE {Herschel OD: 
558). The level 2.5 product for PACS and the level 2 prod¬ 
uct for SPIRE from the pipeline procedures are used for our 
data analysis. Source photometry is carried out using DAOphot 
algorithm in the Herschel Interactive Processing Environment 
(HIPE). We apply beam correction, colour correction, aperture 
correction for a spectral index of -2 and adopt a Hux calibration 
eiTor of 5% at PACS bands and 7% at SPIRE bands as recom¬ 
mended in the PACS and SPIRE Observers Manual. The full 
width at half power (EWHP) beam sizes are 6.8" at 100 yum, 
11.4" at 160 yum, 17.6" at 250 yum, 23.9" at 350 yum and 35.2" 
at 500 yum, respectively. 

3. Results 

3.1. Radio emission from ATCA observations 

In Fig.[TJa) we present the radio continuum emission images at 
20 cm from the ATCA. Among the four LABs, B6 and B7 are 
detected with huxes of 67+17 yuJy and 77+16 /rJy, respectively, 
and B5 is marginally detected at 3 cr (51 + 16 yuJy). Eor all de¬ 
tected sources, their positions are consistent with the central 
positions of the LABs. Only B1 is not detected by the observa¬ 
tions. 
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Fig. 1. ATCA 20 cm, Spitzer MIPS 24jum and Herschel PACS and SPIRE data for the four Lya blobs (LABs) in J2143-4423. a) 
Contours and gray scale maps of ATCA radio emission. The contours are -2, 2, 3, 4, 5 and 6 x 15 piy (1 cr), with a synthesized 
beam of 7.8"x4.8", which is shown in the lower left corner of each panel, b) Gray maps of Spitzer MIPS 24 pm emission 
(Colbert et al. 2006). c-g) Contours and gray scale maps of Herschel FIR emission. The contours are -2cr, 2cr, 3cr, 4cr, 5cr and 
6cr (see § 2.2 for the noise level of each band). A circle with a diameter of 40" is shown in each panel. The circles in B7 are 
on an off-center position (5", 0") to cover most FIR emission. All sources are centered on the positions of the four LABs (see 
Colbert et al. 2006) as shown with plus signs in each panel. All offsets are relative to the positions of the LABs. 


3.2. FIR emission from Herschel observations 

All four LABs are observed with Herschel PACS at 100 and 
160 pm and SPIRE at 250, 350 and 500 pm, and the images 
are shown in Fig. [TJc-g). The observed flux densities are cal¬ 
culated for the areas within the blue circles as shown in Fig. [T] 
and are listed in Table [T] B1 is not detected but contaminated 
by a nearby strong source about 20" in the north-west, which 
is the background QSO LBQS2138-4427 at z = 3.2 (Francis & 
Hewett 1993), and its emission features at different FIR bands 
appear to reach out to B1 from this location. There is no FIR 
counterpart for B5 in any Herschel band. 

3.3. Redshifts of the FIR sources 

To estimate the redshift of the FIR sources associated with B6 
and B7, we try to ht the data with the SEDs of different tem¬ 
plates (Polletta et al. 2007) at different redshifts and hnd that 
the starburst templates can well reproduce the data. With the 
observational data and the SEDs of the templates, the mini¬ 
mum reduced value for each redshift can be calculated and 
the corresponding probability can be estimated. In this analy¬ 


sis, we include hve Herschel band, APEX 870 pm data (Beelen 
et al. 2008), and Spitzer MIPS 24 pm data (Colbert et al. 2006). 

Among four typical templates, Arp 220, M 82, Mrk 231 and 
NGC 6240, we hnd that the spectral energy distribution of star- 
burst galaxies NGC 6240 and Arp 220 ht the data best, and Mrk 
231 doesn’t ht well because it has warm IR emission from its 
AGN which is not really consistent with the data. Fig. shows 
the probability distribution against redshift for both LABs. The 
estimated redshifts are 2.20 ^q3j for B6 and 2.20 ;^q3q for B7 
respectively. Considering the uncertainty of this method to de¬ 
termine the redshifts, both values are consistent with the Lya 
redshift of 2.38 of the LABs. Adopting the number count study 
of Herschel sources in Clements et al. (2010), the probability 
of hnding a 350 pm source with a hux greater than 40 mJy 
within 20 arcsec is 2%. According for such a low number den¬ 
sity of strong FIR sources and the positional coincidence of the 
LABs with strong FIR sources, the FIR sources are very likely 
associated with the LABs. Nevertheless, future spectroscopic 
observations from molecular lines at millimeter or from for¬ 
bidden lines at near-infrared will be quite important to conhrm 
it. In the following sections, the Lya redshift of 2.38 will be 
adopted for the LABs. 
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Table 1. Observational and derived parameters towards the four LABs“ 


Source 

20 cm* 

[Fly] 

100 pm 
[mJy] 

160 pm 
[mJy] 

250 pm 
[mJy] 

350 pm 
[mJy] 

500 pm 
[mJy] 

[lO'^Le] 

Mdust 

[lO'Mo] 

B1 

<51 

<4.2 

<9.0 

<17.9 

<19.6 

<22.5 

<2.8 


B5 

51+16 

<2.1 

<11.1 

<17.5 

<18.7 

<19.8 

<2.5 


B6 

67+17 

13.2+3.2 

53.9+8.0 

49.7+9.0 

53.7+10.7 

36.7+10.3 

10.0+1.9 

3.2+0.6 

B7 

77+16 

12.9+4.0 

33.5+10.0 

41.6+7.8 

48.0+10.6 

39.2+8.6 

8.6±2.3 

5.0+1.0 


“ The wavelengths shown in this table are the redshifted values. 

Measured fluxes have been modified by a primary beam correction (less than 15%). 

The total luminosities are calculated between rest frame wavelengths of 40 pm to 200 pm from the dust models (see § |3.4| for details). 
The 3 cr upper limits are given for undetected sources. 
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Fig. 2. Probability as a function of redshift for B6 and B7. NGC 6240 and Arp 220 are adopted as the most appropriate starburst 
templates for B6 and B7, respectively. A red vertical line denotes a redshift of 2.38 for Lya emission. 


3.4. Dust properties 


For B6 and B7, we have included the measurements from the 
hve Herschel bands as well as the 870 pm data taken from 
Beelen et al. (2008) in the dust continuum analysis using 
a single-component dust model as described in WeiB et al. 
(2007). Spitzer MIPS 24 pm data (Colbert et al. 2006) are not 
used in the model htting because they are strongly affected by 
PAH features, but are shown in Fig. to allow for a better 
comparison with overlaid templates. We hnd a dust tempera¬ 
ture, Tdust, of 70+5 K and a dust mass, Mdust, of (3.2+0.8)xl0^ 
Mo for B6, and Tdust = 70+5 K and Mjust = (5.0+1.0)xl0** Mq 
for B7, respectively. The implied FIR luminosities are LpiR 
= (10.0+1.9)x10'2Lo forB6, and Lfir = (8.6+2.3)xl0'2 Lg 
for B7, respectively, where Lrir is integrated from 40 pm to 
200 pm in the rest frame. The upper Lrir limits for both B1 
and B5 are ~2.5-2.8x10'^ Lo- 


3.5. Star formation rates 

Here we derive the star formation rates from the Lya, far- 
infrared and radio luminosities. To estimate the star forma¬ 
tion rate (SFR) from the Lya luminosity, we first assume that 
star formation (SF) powers the observed Lya flux. We use 
an unreddened Lya/Ha ratio of 8;1 and the conversion fac¬ 
tor between Ha luminosity and SFR (Kennicutt 1998), yield¬ 
ing SFR(Lya)/(Mo/yr) = LLyQ’/(10"^^ erg s“'). This provides a 
lower limit because the extinction of Lya emission caused by 
dust will largely reduce the observed Lya luminosity. With the 
FIR luminosity derived from Herschel data, we can estimate the 
SFR by using the relation SFR(LFiR)/(Mo/yr) = 1.7 xLfir/(10*° 
Lq) (Kennicutt 1998). If the observed radio emission, with 
a rest wavelength of 6 cm, is dominated by free-free emis¬ 
sion in HII regions, one can also relate the SFR by the re¬ 
lation SFR(Li.4GHz)/(Mo/yr)- 5.52x10-22 Li.4ghz/(W Hz-') 

(Bell 2003). The radio luminosity at 1.4 GHz at the rest frame 
can be estimated from the observed flux at 1.51 GHz by assum¬ 
ing a relation S oc v“, where S is the flux density and the typical 
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Fig. 3. Single-component dust models for B6 and B7 (a redshift of 2.38 is adopted). The black solid lines show the thermal 
dust continuum emission of the 70 K dust components for both B6and B7. The open circles represent the measurements at five 
Herschel bands in this paper and the filled circles indicate the flux densities at 24 fim (Colbert et al. 2006). The filled square 
denotes the flux density (or its upper limit) at 870 ^m, taken from Beelen et al. (2008). The wavelengths at the rest frame are 
labelled on the top. For the single-component dust model s ado pted in the figure (see § 3.4 for details of the dust models.), the 
values are 1.1 for B6 and 1.0 for B7, respectively. In § 


3.3 


four typical starburst templates, NGC 6240, M 82, Mrk 231 and 
Arp 220 (Polletta et al. 2007), are adopted to estimate the redshifts for B6 and B7, and their best fits are overlaid in colored lines. 


spectral index a of -0.8 is commonly adopted for the SMGs 
(e.g., Ivison et al. 2010). These values are listed in Table|^ 


Table 2. Derived star formation rates towards the four LABs 


Source 

SFR(Lfir) 

[Mo/yr] 

SFR(L,.4GHz) 

[Mo/yr] 

log LLya" 
[ergs s^‘] 

SFRlLyff) 

[Mo/yr] 

B1 

<480 

<1090 

43.9 

79 

B5 

<430 

1090±340 

43.8 

63 

B6 

1700±320 

1430±360 

43.8 

63 

B7 

1460±390 

1650±340 

43.5 

32 


“ The Lya luminosities are adopted from Colbert et al. (2006). 


4. Discussion and Conclusions 

A high detection rate of radio emission (three out of four) 
around LABs suggests that most LABs do not originate from 
cooling radiation. Instead, photoionization from starbursts 
and/or AGNs may power the LABs in most cases. The high rate 
of FIR detections (two out of four) points to a star-formation 
origin of the LABs. The SEDs of B6 and B7 can be well de¬ 
scribed by starburst dominated templates, as shown in Fig. 


further supporting Lya emission related to the SF in the LABs. 
In B6 and B7, the SFRs derived from Lya fluxes are far be¬ 
low those estimated from FIR luminosities (Table|^. This sug¬ 
gests that the dust indeed greatly reduces the measured Lya 
flux. Comparing the different SFRs, the dust absorption opti¬ 
cal depth of the Lya emission becomes ~3.1-3.6. The SFRs 
estimated from the FIR and radio luminosities are compara¬ 
ble, indicating that the radio emission is dominated by SF, not 
by AGNs. The energetic starbursts can provide enough ioniz¬ 
ing photons to ionize neutral hydrogen atoms in the interstel¬ 
lar medium (ISM), and each subsequent recombination has a 
probability of ~ 2/3 of ending up as a Lya photon (Partridge 
& Peebles 1967). After escaping the galaxy’s ISM, these Lya 
photons can be resonantly scattered by neutral hydrogen atoms 
in the intergalactic medium (IGM), which tends to make the 
Lya emission extended (Zheng et al. 2011). 

Cen & Zheng (2013) propose an SF-based model and pre¬ 
dict that LABs at high redshift correspond to protoclusters con¬ 
taining the most massive galaxies and cluster halos in the early 
universe as well as ubiquitous strong infrared sources under¬ 
going extreme starbursts. This may be supported by the multi¬ 
ple Spitzer/MIPS sources detected in both LABs (see Fig[2b), 
Colbert et al. 2006, 2011). Indeed, Prescott et al. (2012b) sug¬ 
gest that LABs may be the seeds of galaxy clusters by resolving 
















6 


Y. Ao et al.: What Powers the Lya Blobs? 


the galaxies within a LAB at z = 2.7. The strong FIR emission 
and the inferred high SFRs support the presence of a strong 
starburst in both B6 and B7. However, AGN-dominated tem¬ 
plates like Mrk 231 can not well reproduce the data (see § |3.3| l, 
suggesting that the SF instead of AGN may power the Lya 
emission in both LABs. The model also predicts that the most 
luminous FIR source in each LAB is likely representing the 
gravitational center of the protocluster. Fig. [TJc-g) shows that 
the FIR emission indeed peaks in the centers of B6 and B7. The 
radio continuum emission is detected exclusively in the centers, 
which suggests that the source with most luminous FIR emis¬ 
sion (therefore highest SFR) is in the gravitational center of 
each LAB. Another very important prediction of this model is 
that the Lya emission from photons that escape the galaxy are 
expected to be significantly polarized, which has been for the 
first time confirmed towards LABI in the SSA22 field by Hayes 
et al. (2011), supporting models with central power sources. 
Adopting a gas-to-dust mass ratio of 150 and the SFRs esti¬ 
mated above, the timescales of B6 and B7 are relatively short 
(~100 Myr), which is much shorter than the galaxy building 
timescale. Note that this timescale is a lower limit because (1) 
the LABs may have been alive for a while now, and (2) addi¬ 
tional gas may be continuously accreted. In any case, the LABs 
are visible only during a short time interval during the lifetime 
of their parent clusters. 

Note that the so-called “SF-based model” proposed by Cen 
& Zheng (2013) also includes AGN powering or any central 
powering. The morphologies of the Lya emission of the four 
LABs are quite different (Palunas et al. 2004); B1 and B5 have 
core-like structures, while B6 and B7 are characterized by dif¬ 
fuse and extended emission with physical sizes of ~60-70 kpc. 
The latter may be driven by multiple sources as suggested by 
the MIPS data and are consistent with the SF-based model. 
There is no clear FIR emission detected around B1 and B5. 
Therefore, the Lya emission in both LABs is unlikely predom¬ 
inantly triggered by SF. Overzier et al. (2013) conclude that 
in B1 the photoionization from an AGN is the main driver of 
Lya emission. However, Francis et al. (2013) shows that the 
observed Lya emission in B1 is of complex origin, dominated 
by the sum of the emission from the sub-haloes where the cold 
gas is being lit up most likely by a combination of tidally trig¬ 
gered star formation, bow shocks, resonant scattering of Lya 
from the filament collisions and tidal stripping of the gas. In B5 
radio emission is tentatively detected and therefore the AGN 
may also power the Lya emission. Among the four LABs in 
J2143-4423, two of them, B6 and B7, are mainly driven by 
SF. However, the other two LABs, B1 and B5, without clear 
FIR detection, are predominantly driven by the AGNs or other 
sources of energy still to be specified, but not mainly by star 
formation. We thus conclude that LABs must be powered by 
quite diverse sources of energy. 

With its high angular resolution and superb sensitivity, fu¬ 
ture observations with the Large Atacama Millimeter Array 
(ALMA) will reveal more details about the nature of LABs 
such as testing the predictions of models where the ionization 
is provided by intense star formation and confirming the signif¬ 
icantly polarized dust emission at mm/submm wavelength. 
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